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ABSTRACT 

We present non-linear hydrodynamic pulsation models for OGLE-BLG-RRLYR-02792 
- a 0.26 Mq pulsator, component of the eclipsing binary system, analysed recently by 
Pietrzynski et al. The star's light and radial velocity curves mimic that of classical 
RR Lyrae stars, except for the bump in the middle of the ascending branch of the 
radial velocity curve. We show that the bump is caused by the 2:1 resonance between 
the fundamental mode and the second overtone the same mechanism that causes the 
Hertzsprung bump progression in classical Cepheids. The models allow to constrain 
the parameters of the star, in particular to estimate its absolute luminosity (« 33 L ) 
and effective temperature 6970 K, close to the blue edge of the instability strip). 

We conduct a model survey for the new class of low mass pulsators similar to 
OGLE-BLG-RRLYR-02792 - products of evolution in the binary systems. We compute 
a grid of models with masses corresponding to half (and less) of the typical mass 
of RR Lyrae variable, 0.20 M Q ^ M < 0.30 M©, and discuss the properties of the 
resulting light and radial velocity curves. Resonant bump progression is clear and 
may be used to distinguish such stars from classical RR Lyrae stars. We present the 
Fourier decomposition parameters for the modelled light and radial velocity curves. 
The expected values of the ip^i Fourier phase for the light curves differ significantly 
from that observed in RR Lyrae stars, which is another discriminant of the new class. 

Key words: hydrodynamics - methods: numerical - binaries: eclipsing - stars: os- 
cillations - stars: variables: RR Lyrae 



1 INTRODUCTION 

OGLE-BLG-RRLYR-02792 (below shortened to RR) is a 
large amplitude, single-periodic (Pq « 0.6275 d) pulsator, 
component of a physical, well detached, double-lined eclips- 
ing binary system discovered in the Galactic bulge by the 
Optical Gravitational Lensing Experiment (OGLE) survey 
l|Soszvnski et al.ll201ll ) and analysed in detail i n the course 
of the Araucaria project by IPietrzvriski et al.l ([2012) . Dis- 



entangled /-band light curve and radial velocity curve are 
presented in Fig. [T] The pulsation period and shape of the 
light curve places the star among classical RR Lyrae stars 
pulsating in the fundamental mode (RRab stars), which is 
clear from the analysis of the Fourier decomposition param- 
eters, amplitude ratios, Rki, and Fourier phases, tpui- They 
are constructed through the fit of the Fourier series to the 
data, i.e., 
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mi = A + A k sin(fejt + tp k ), 
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and then, 

Rki = Ak/Ai, ipki = <pk — k<pi, (2) 

where ui = 2tt/P is pu lsation frequency. As computed by 
iPietrzvriski et ail (|2012t ) the pulsation period decreases with 
the rate of (8.4±2.6) x 10 -6 dyr -1 . Here we adopt the mean 
pulsation period derived from the /-band photometry, which 
spans nearly 10 years, Po = 0.627527 d. Fourier decomposi- 
tion parameters of RRLYR-02792 are given in Tab. Q] and 
plotted in Fig. [2] along with the OGL E data for thousands o f 
RRab stars from the Galactic bulge (|Soszvriski et al.ll201ll ). 
Note that the Fourier phases from the OGLE catalogue, 
given in the cosine convention, were shifted to match the 
sine convention (eq. 1) that we consistently use in this pa- 
per (for both the light and radial velocity curves Fourier 
phases of RRLYR-02792 fit the typical values of RR Lyrae 
stars very well. Amplitude ratios are low, but fit the tail well 
visible in the progression for RRab stars between Po = 0.6 
and 0.65 days. 

Surprising ly, the dynamical mass of RRLYR-02792 de- 
rived by IPietrzvriski et alj (|2012T ) is very low, M = 0.261 ± 
0.015 M©, which is roughly two times smaller than the typ- 
ical mass of the RR Lyrae star. It is too low to initi a te the 
helium burning. As discussed by IPietrzvriski et al.l (|2012r i 
the star is a result of mass transfer in the binary sys- 
tem. Other parameter s of RRLYR-02792 relev ant for this 
study and derived by IPietrzvriski et al] (|2012T l are radius, 
R = 4.24 ± 0.24 R© and estimate of the effective tempera- 
ture, T e ff = 7320 ± 160 K (rough estimate based on the as- 
sumed 5000 K temperature of the secondary 1.67 Mq giant 
component and temperature ratio derived from the analy- 
sis of the light curve). Corresponding absolute luminosity 
may be estimated to be L = 46.2 ± 9.3 Lq. No metallicity 
determination is available for RRLYR-02792. 

RRLYR-02792 is a new member of class of pulsators, 
for which mass transfer associated with binary evolution has 
significant effect on their pulsation properties. We will call 
such stars binary evolution pulsators (BEPs in the follow- 
ing) and we will call the RRLYR-02792 the BEP. Other well 
known examples of the class are hot su bdwarf pulsators , 
members of eclipsing binary systems fsee iFor et al.ll2~010h . 
We also note that BEP may be a single object that went 
through binary evolution in the past, and was ejected from 
the binary because of e.g. tidal interactions in the globular 
cluster. 

In this paper we present the first pulsation calcula- 
tions for the BEP, both linear and non-linear, full- amplitude 
models. Nearly exact fit of the light and radial velocity 
curves is possible, which allows to constrain the stellar pa- 
rameters and provides an insight into the pulsation dynam- 
ics of the star (Section|3}. We also present a non-linear model 
survey with parameters of the models corresponding to the 
expected parameters of the members of the new class, i.e. 
with masses 2 — 3 times smaller than typically assumed for 
RR Lyrae stars (Section [4j . The main goal is to provide the 
criteria to distinguish the BEPs from classical RR Lyrae 
stars, for which binary nature cannot be inferred directly 
because of too low inclination angle. In this case, we only 
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Table 1. Fourier decomposition parameters for the /-band light 
curve and radial velocity (RV) curve of RRLYR-02792 (the BEP, 
sine decomposition). Standard errors are given in the parentheses. 
In the last line the peak-to-peak amplitude, A, is given. 
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Figure 1. Disentangled /-band light curve (top) and radial ve- 
locity curve (bottom) for RRLYR-02792 (the BEP). Solid lines 
represent the 10th (light curve) and 6th (radial velocity curve) 
order Fourier fits to the data. Fourier decomposition parameters 
are given in Tab.fT] Note that we use an observers' convention in 
the plot of radial velocity curve, i.e. negative values correspond 
to expansion. 

see the light/radial velocity curve of a pulsating component, 
which, as in the case of the BEP, may strongly mimic the 
typical curve of RR Lyrae pulsator. 



2 NUMERICAL METHODS 

All our models are computed with the convective pulsa- 
tion codes described in detail bv lSmolec fc Moskalikl |2008). 
These are static model builder, linear non-adiabatic code 
and non-linear direct time-integration hydrodynamic code. 
The codes are Lagrangian. Radiative transfer is treated with 
a simple diffusion approximation. For the convective energy 
transport we use the tim e-dependent, turbulent convection 
model of iKuhfufj I 19861 ). It is a simple one-dimensional, 
one-equation formulation that contains several dimension- 
less scaling parameters. These are the mixing- length pa- 
rameter, a, and parameters scaling the turbulent fluxes and 
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Figure 2. Pcak-to-peak amplitude, A, and Fourier decomposi- 
tion parameters for the //-band light curves of the Galactic bulge 
RR Lyrae stars of the OGLE survey (grey dots, each third plotted; 
ISoszvnski et alj|201lh . RRLYR-02792 (the BEP, black diamond) 
and hydrodynamic models adopting convcctive parameters of set 
SI (lines, see the keys in the top and bottom panels). 

the terms that drive/damp the turbulent energy: a p (turbu- 
lent pressure), a m (eddy- viscous dissipation), q c (convective 
heat flux), at (kinetic turbulent energy flux), a s (buoyant 
driving), ad (turbulent dissipation) and 7 r (radiative cool- 
ing). Theory provides no guidance for their values. For a B , 
a c , «d and 7 r a values which result from comparison of the 
static time-independent version of the model with the Mix- 
ing Length Theory (MLT) are in use (see caption of Tab. [2]). 
In practice, the values of convective parameters are cali- 
brated to match as many observational constraints as pos- 
sible. Slightly different values are necessary for modelling 
different groups of stars, e.g. Cepheids and RR Lyrae stars 
of different stellar systems. It makes the modelling of BEPs 



particularly difficult. The only known member of the class, 
the BEP, provides very little constraints on the models. 
Consequently, several sets of values of convective parameters 
need to be considered - Tab. [2] As a first guess we adopt the 
parameters that well reproduce the observed properties of 
the Galactic bulge RR Lyrae stars, in particular the Fourier 
parameters of the /-band light curves. In Fig. [2] we com- 
pare the computed Fourier parameters for models adopting 
convective parameters of set SI, and assuming typical pa- 
rameters of RR Lyrae stars (masses 0.55 Mq and 0.65 Mq, 
luminosities: 40 — 70 Lq, metallicities, Z: 0.01, 0.001, and 
0.0001) with OGLE observations. We get a good agreement 
for amplitudes, Fourier phases and R21 ratio and somewhat 
too low model values for R31. Similar agreement (with much 
better match of R31) is obtained for models of set S2, which 
include the effects of radiative cooling. In set S3 the con- 
vective parameters are the same as in set SI, except eddy 
viscosity which is lower. Consequently, we expect higher pul- 
sation amplitudes in models adopting parameters of set S3. 
In linear model survey we also consider set S4, in which 
effects of turbulent pressure are turned on. Inclusion of tur- 
bulent pressure in non-linear models leads to long-lasting in- 
tegration and convergence difficulties for some models. Con- 
sequently, in this first non-linear model survey for BEPs, we 
neglect the effects of turbulent pressure and consider models 
of sets SI, S2 and S3, only. 

Our models are chemically homogeneous, non-rotating, 
envelope models. They consist of 150 mass shells extend- 
ing down to 2 x 10 6 K. It is a standard depth in modelling 
the Cepheid and RR Lyrae pulsation. It is deep enough 
to cover the region relevant for k mechanism driven pulsa- 
tion but too shallow to enter the inhomogeneous hydrogen 
burning shell region. In order to capture and resolve the 
thin hydrogen ionisation region the anchor zone of temper- 
ature equal to 11 000 K is located 40 zones below the sur- 
face. Non-linear model integration is conducted till the limit 
cycle pulsation (single periodic, full-amplitude pulsation) is 
reached, which requires computation of 1 000 up to 5 000 pul- 
sation cycles. Each pulsation cycle is covered with 600 time- 
steps. The bolometric lig ht curve s are t ransformed to the 
/-band through the static iKuruczl (|2005l ) atmosphere mod- 
elfl Projection factor for the BEP is unknown. We scaled 
the model velocity curves by projection factor used fo r clas- 
sical RR Lyrae pulsators, i.e. p = 1.38 l|Fernlevlll994l ). 

The chemical structure of the star is subject of uncer- 
tainty. No metallicity determination is available. Also, no 
realistic evolutionary computations exist for this type of bi- 
nary evolution. Significant fraction of the envelope was re- 
moved during the evolution. Pulsation is driven in the out- 
ermost layers which, although very extended, contain only a 
fraction of percent of the total stellar mass (typically be- 
low 1 per cent in our models) . In this study we assume 
that the outer shell has preserved the standard chemical 
composition typical for Population II stars, i.e., X = 0.76, 
Z = 0.01 — 0.0001. We also compute model sequences with 
X — 0.70 which corresponds to higher helium content as 
in Popula tion I stars. In all our models we adopt the OP 
opacit ies (|Seatonll2005h and solar mixture of lAsplund et all 
|2004h . 
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Figure 3. The HR diagrams for 0.26 Mq models adopting different values of convective parameters and different chemical compositions 
as indicated in the top right part of each panel. Instability strips are indicated with gray-shaded areas (light-gray for the fundamental 
mode, dark-gray for the first overtone). Solid line is the line of constant fundamental mode period, Pq = 0.6275 d. Long-dashed lines 
indicate the radius determination for the BEP (with la errors) and dotted line shows the loci of the 2:1 resonance with the second 
overtone, P2/P0 = 0.5. Circle marks location of the static model from Fig. [4] and diamonds mark the location of the best matching 
models for the BEP (Section [372j Tab. |3j. 



Table 2. Convective parameters adopted in the models. Parame- 
ters a B , ct c , oy, dp and 7 r are given in the units of standard values 
(q s = q c = 1/2^/273, a d = 8/3^273, a p = 2/3 and -y r = 2V3; 
see lSmolec fc Moskalik 2008|, for details). Only first three sets of 
parameters were used in non-linear computations. 
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3 PULSATION MODELS FOR THE BEP 

In this Section we focus on modelling the observed properties 
of the BEP. We first conduct a linear model survey and 
discuss properties of the static models (Section 13. ip . Next 
we focus on non-linear modelling of the /-band light curve 
and radial velocity curve of the BEP (Section l3.2p . 



3.1 Linear model survey 

The models are computed along horizontal stripes of con- 
stant luminosity (L — 16 — 40 Lq, step 2Lq) and with 100 K 



step in effective temperature. We adopt M = 0.26 Mq for all 
the models. Computed Hertzsprung-Russel (HR) diagrams 
for different sets of convective parameters and assuming dif- 
ferent chemical compositions are plotted in Fig. [3] Except 
for model computations adopting convective parameters of 
set S2, first overtone is linearly damped in our models. The 
typical width of the fundamental mode instability strip is 
1 000 K and increases towards the higher luminosities. 

Before we discuss the models appropriate for the BEP, 
we first discuss general properties of typical static model, 
as pulsation properties of BEPs were not studied so far. In 
Fig. |3] we plot results for a model located in the middle of 
the instability strip (T e g = 6500 K, L = 25 Lq, marked with 
open circle in the top left panel of Fig. [3j . Its period is only 
slightly longer than the period of the BEP. The top panel 
shows the convective heat flux (solid line) and adiabatic tem- 
perature gradient (dashed line). In the middle panel of Fig. [4] 
we plot the logarithm of Rosseland opacity (solid line) and 
its logarithmic temperature derivative (dashed line) - cru- 
cial quantities for studying the pulsation driving. The metal, 
helium and hydrogen opacity bumps are marked with labels. 
In the bottom panel the work integrands are plotted - total 
with the solid line and contribution of eddy viscous work 
with the dashed line. Pulsation driving occurs in regions 
with positive work integrand. Eddy viscosity always damps 
the pulsation and its contribution is negative. Model proper- 
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ties are very similar to the pro perties of classical RR Lyrae 
models (e.g. IStellingwerj|l982l V It is clear that pulsation is 
driven through the k mechanism working in the hydrogen- 
helium ionisation zone. We note a significant contribution 
from the hydrogen ionisation region. The main envelope con- 
vection zone is associated with the hydrogen-helium (Hel-H) 
ionisation region. The second, less effective convective region 
develops deeper at the Hell ionisation zone. All ionisation 
regions are properly resolved in our models as run of V a d 
indicates. 

In each panel of Fig. [3] we plot a line of constant funda- 
mental mode period, Po = 0.627527 d. The dashed lines in 
Fig.[3]are lines of constant photospheric radius R — 4.24 Rq, 
with la errors (±0.24 Rq), which correspond to the values 
measured for the BEP. For models with mass corresponding 
to the measured mass of the BEP, these two lines (of con- 
stant Po = 0.6275 d and of constant R = 4.24 Rq) should 
roughly overlap (period-mean density relation), which is not 
the case for any of the considered model parameters. The 
lines are nearly parallel but shifted with respect to each 
other. The best agreement is obtained for models of the 
lowest metal abundance and for the models including the 
effects of turbulent pressure, which expands the model en- 
velope (middle bottom panel in Fig. [3} • In this case the dif- 
ference is la (at high L). For other models the difference is 
larger, however, always smaller than 2a. Hence, we treat this 
discrepancy as a warning sign and note that multi-band pho- 
tometric observations are necessary to derive the radius of 
the BEP more accurately (see also discussion in Section [5}. 

In the HR diagrams in Fig. [3] we also plot the loci 
of the 2:1 resonance between the fundamental mode and 
the second overtone, P2/P0 = 0.5 (dotted lines). This 
resonance causes the famous Hertzsprung bump progres- 
sion in classi cal Cepheids pulsating in the fundamental 
mode (e.g. Simon fc Schmidt]! 19761 ; iKovacs fc Buchlerll 19891 ; 
iBuchler. Moskalik fc Kovacsl 11990), In the discussed mod- 
els it lays well within the instability strip, independent of 
the model parameters (see also Fig. [5] for models with dif- 
ferent mass). Our non-linear computations presented in the 
next sections show that this resonance strongly influences 
the pulsation of BEPs. 



3.2 Non-linear model survey 

The non-linear models for the BEP are computed along lines 
of constant fundamental mode period, Po = 0.6275 d. Mod- 
els are labelled by their absolute luminosity, which increases 
in a step of ILq. All models have M = 0.26 Mq. The least 
luminous models are located at the red edge of the instabil- 
ity domain and, as luminosity increases, shift towards the 
blue (see Fig. [3} . We computed six model sequences: with 
convective parameters of sets SI, S2 and S3 and with metal- 
licities Z = 0.01 and Z = 0.001 {X = 0.76) in each case. For 
models with convective parameters of set S3 we also com- 
puted two model sequences with X = 0.70 (Z = 0.01 and 
Z = 0.001). In Fig. [5] we show a full set of the computed I- 
band light curves and radial velocity curves for models with 
convective parameters of set SI and with Z = 0.01. At phase 
the radius reaches its maximum value (velocity is equal to 
0). The radial velocity curves have a usual triangular shape. 
One may note a systematic progression of the curve's shape 
as luminosity changes. For the lowest luminosity models the 
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Figure 4. Properties of the exemplary static model with M = 
0.26 Mq , L = 25 L Q and T cff = 6500 K (marked with circle in the 
top left panel of Fig. [3J). Top panel: convective heat flux (solid 
line) and adiabatic temperature gradient (dashed line); middle 
panel: logarithm of opacity (solid line) and its logarithmic tem- 
perature derivative (dashed line); bottom panel: work integrands, 
total (solid line) and eddy- viscous contribution (dashed line). 



velocity curves are more or less symmetrical, the descend- 
ing branch is even slightly longer than the ascending branch 
for the model with L = 19 Lq. The velocity varies slowly at 
the phase of maximum expansion velocities. For models with 
L = 22 — 23 Lq the velocity curve is nearly flat at this phase. 
As luminosity increases a bump appears on the ascending 
branch (visible already for L = 25 Lq at phase around 0.8), 
which moves upwards towards the phase of maximum radius 
(L = 31 Lq) and then, vanish (L = 32 Lq). 

Considering the /-band light curves, at the lowest lumi- 
nosities they differ significantly from the typical RR Lyrae 
curves. A pronounced bump is present on the ascending 
branch - two brightness maxima are observed during one 
pulsation cycle. As luminosity increases the bump becomes 
less pronounced, turns into narrow standstill (28, 29 Lq) 
and vanishes. At highest luminosities, the light curves re- 
semble those of typical RR Lyrae stars, with the except for 
less sharp light variation in the brightness maximum. 

Described variation of the shape of radial velocity and 
/-band light curves suggests, that the 2:1 resonance with 
the second overtone is influential here, and shapes the bump 
progression. It is not very pronounced because we consider 
a specific model sequences of constant period. It is clear 
form Fig. [3] that most of our models lay to the blue of 
the P2/P0 = 0.5 line and thus, majority of models have 
P2/P0 > 0.5 - see exact values on the right side of Fig. [5] 
In the case of bump progression observed e.g. in classical 
Cepheids, which is caused by the same resonance, the shape 
of the light /radial velocity curves is affected by the reso- 
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Figure 5. /-band light curves (left panel) and radial velocity curves (right panel, observers convention, i.e. negative values correspond 
to envelope expansion) for hydrodynamic models of 0.26 Mq computed along a line of constant fundamental mode period corresponding 
to the pulsation period of the BEP. Models adopt convective parameters of set SI. For clarity, the consecutive light curves were shifted 
by 0.12 mag and radial velocity curves by lOkms -1 . Velocity curves were scaled by constant projection factor, p = 1.38. Thick red lines 
indicate the best matching model for the BEP. 



nance in a wide range of period ratios, 0.48 < P2/P0 < 0.54. 
The position of the bump depends critic ally on the value of 
P2/P0, as model computations indicate (iKovacs fe Buchlerl 
1 19891 : iBuchler. Moskalik fc Kovacsl[l99oh . For P2/P0 > 0.5 
the bump is located on the ascending branch of the velocity 
curve, which is also the case in our models. The full bump 
progression in our models will be demonstrated in Section [4] 
where we consider models covering the full instability strip 
and thus probing the P2/P0 ratio in a wide range (Fig. 1 11 p . 

It is clear that among the models displayed in Fig. [S] 
these are the highest luminosity ones that reproduce the 
observations of the BEP best (cf. Fig. [TJ. in particular, the 
model for L = 31 L© (marked with thick red lines in Fig. [5} 
seems to be the best one. A more quantitative selection of 
the best models for the BEP is based on comparison of the 
observed and computed Fourier decomposition parameters 
of the lowest order, the amplitude Ai, amplitude ratios, R21 
and R31, and Fourier phases, Lp2i and (f3i, for both /-band 
and radial velocity curves. In Figs. [6] and [7] we show such 
comparisons for specific Fourier parameter planes, i.e., (f2i — 
R21 and (£31 — R31 and three sequences of models. We stress 
that periods of all the models match the period of the BEP. 
Models differ in luminosity (and effective temperature) - 
the least luminous model (and the coolest one) is marked 
with open circle. Luminosity of the consecutive models in 
a sequence differs by 1 Lq . Again it is clear that these are 



the highest luminosity models that match the observations 
of the BEP best. 

Yet another observable that can be compared with the 
models is the phase lag between the radial velocit y curve 
and t he light curve, defined as A<pi = (p\ al — (p™ as (jSimonl 
1 1984! ) . For the BEP we get A<^i = 0.055 ± 0.02/0. Unfor- 
tunately there are no phase lag determinations for classical 
RR Lyrae stars using the /-band light curves. With the V- 
band light curves, the mean p hase lag for a sample of RR 
Lyra e stars is {Atpi) = —0.19 ()Ogloza. Moskalik fc Kanburl 
2000) . It is negative for all the star s used in the work 
of lOgloza. Moskalik fc Kanburl (|2000l ) (Moskalik private 
comm.). Our models for classical RRab stars (displayed in 
Fig. [2]) point however, that phase lags computed using the 
/-band light curves are systematically larger than phase lags 
computed using the V-band light curves and may have both 
positive and negative values. 

We now compute the offsets between the model and the 
observations in the Fourier parameter space, defined as 



3 Because of the pulsation period change, A<pi was computed 
using the last 500 days of the OGLE-III photometry only. At the 
end of this span the radial velocities were measured. The adopted 
pulsation period (P = 0.627510(5) d) was derived from the pho- 
tometry alone, but well agrees with the (less accurate) period 
derived from radial velocity measurements. 
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Figure 6. Plots of the Fourier parameters, ip2l vs. R21 (top) 
and </?3i vs. R31 (bottom) for the computed /-band light curves 
for three sequences of models of constant period (all models with 
X = 0.76 and Z = 0.01, but with different convective parame- 
ters, as indicated in the key). Open circle marks the first, lowest 
luminosity model in each sequence. Consecutive models differ by 
1 Lq . Larger symbols show the best matching models for the BEP 
(marked with diamond). 



d = ^/^(pi,mod - Pi,obs) 2 /P?,obs> (3) 

where pi denotes one of the low order Fourier parameters 
and amplitude, pi 6 {^4, R21, R31, V21, <P3i}, for our models 
(pi.mod) and for the observed curves (pi,obs, Tab.Q]). Depend- 
ing on the Fourier parameters used, several offsets may be 
constructed. At first we consider three definitions of d. The 
first offset, (fe, is computed in the 2D-plane of the lowest or- 
der Fourier coefficients, </?2i — R21 (top panels of Figs. [6] and 
[7|. The second offset, d,4, is constructed of all four low order 
amplitude ratios and Fourier phases, and the third offset, 
ds, includes also the amplitude. 

We note that for a given model sequence and for both 
the /-band and radial velocity curves the three offsets are 
minimised either for the same model or by models that dif- 
fer by ILq. For example, for models adopting convective 
parameters of set S2, and Z — 0.01 all offsets for the I- 
band light curves and radial velocity curves are minimised 
in the model with L — 32 Lq, except d± for the radial veloc- 
ity curves which is minimised in the model with L = 31 Lq. 
In the more complex case of models with convective param- 
eters of set SI (and Z — 0.01) di and ds for the /-band 
light curves and d2 for the radial velocity curves are min- 
imised in the model with L = 31 Lq, while d2 for the /-band 
light curves and di and ds for the radial velocity curves are 
minimised in the model with L — 30 Lq. 

For the final selection of the best model we use the off- 
set d a u constructed from all low order Fourier decomposition 
parameters and amplitudes for both the /-band light curve 
and radial velocity curve (ten parameters). The model that 



Figure 7. The same as Fig. [6] but for the radial velocity curves. 



minimises d a ii for a given model sequence is considered the 
best one. Properties of our best matching models for all con- 
sidered model sequences are given in Tab. [3] In the last row 
we provide the value of dan. In the HR diagrams in Fig. [3] 
location of the best models is indicated with diamonds (in 
relevant panels). In Figs.[6]and[7]the best models are marked 
with larger symbols. In Fig. |H]the light and radial velocity 
curves for our best models for the BEP are compared with 
observations. The light curves were shifted vertically so the 
mean brightness is zero. All the curves were phased in the 
same way - phase corresponds to the maximum radii (ra- 
dial velocity goes through 0). No other shifts were applied 
to the curves (except the vertical shift between consecutive 
models applied for clarity). 

Analysis of Tab. |3] and Figs. [3JH El and El allows to draw 
some interesting conclusions. 

(i) All the best fitting models are located close to the 
blue edge of the instability strip (Fig. [3J • The effective tem- 
peratures and absolute luminosities lie in a narrow ranges of 
~ 6850—7050 K and 31—36 Lq, respectively. More metal rich 
models (Z = 0.01) are slightly (50 — 100 K) cooler and less 
luminous (2 — 3Lq) than metal poor models (Z — 0.001). 

(it) Radius depends on metallicity and is equal to ~ 
3.88 R for Z = 0.01 and « 3.95 R for Z = 0.001. 

(m) Considering the /-band light curves (left panel of 
Fig-Oi clearly the best fit to the observed light curve is ob- 
tained for models with convective parameters of set S3, and 
with X = 0.70 (for both metallicities, bottom two curves in 
Fig. |HJ). These are the models with the lowest value of d a n. 
A slightly worse fit is obtained for models with X = 0.76. 
Among these models the best ones are those with convective 
parameters of sets S2 and S3 and with the higher metallic- 
ity (Z — 0.01, third and fifth curve from top in Fig. [8]). For 
all these models the radial velocity curves well reproduce 
the observed curve, in particular, a bump on the ascending 
branch is present, although less pronounced than in obser- 
vations. We conclude, that chemical composition typical for 
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Figure 8. /-band (left panel) and radial velocity (right panel) curves for hydrodynamic models reproducing the observations of the BEP 
(black solid line, Fourier fit from Fig. [l} best. Models adopting convective parameters of sets SI, S2 and S3 are plotted with red, green 
and blue lines, respectively. 



Table 3. Properties of the best hydrodynamic models for the BEP. 



set 




SI 




S2 




S3 






X/Z 


0.76/0.01 


0.76/0.001 


0.76/0.01 


0.76/0.001 


0.76/0.01 


0.76/0.001 


0.70/0.01 


0.70/0.001 


L/Lq 


31 


33 


32 


35 


32 


34 


33 


36 


T eff [K] 


6864.2 


6913.8 


6910.5 


7003.0 


6911.4 


6957.3 


6966.7 


7053.6 


R/ Rq 


3.88 


3.94 


3.88 


3.95 


3.88 


3.95 


3.88 


3.96 


P2/P0 


0.530 


0.537 


0.536 


0.545 


0.532 


0.539 


0.532 


0.542 




-0.14 


-0.11 


-0.16 


-0.14 


-0.13 


-0.03 


-0.07 


-0.06 


dull 


0.597 


0.421 


0.382 


0.742 


0.290 


0.630 


0.200 


0.227 



population I stars (X — 0.70) seems more appropriate for 
the BEP (see discussion in Section [SJ. 

(iv) For the lower metallicity models (Z = 0.001) a 
phase shift between the model and the observed light curves 
is apparent for models of all three sets of convective pa- 
rameters (left panel in Fig. [SJ. A minimum brightness in 
these models follows the observed minimum by ~ 0.1 — 0.15 
pulsation cycle. The shift is also present in the case of the 
radial velocity curves, for all the models. A phase of max- 
imum expansion velocity in our models occurs later than 
is observed. Again the difference is slightly larger for the 
lower metallicity models. To explain these phase shifts we 
first note that we defined phase as the phase of maximum 
radius. It is determined with the help of the radial veloc- 
ity, which changes the sign. At around this phase a bump 
is present in the radial velocity curves, which manifests it- 
self as a slower (less steep) variation of the radial velocity. 
Thus, its exact location may influence our phasing of the 
model curves. The location of the bump depends on the 
P2/P0 period ratio which is a function of metallicity, as is 



evident from an analysis of the values given in Tab. [3] At 
lower metallicity, the P2/P0 ratio becomes larger and the 
bump moves up the ascending branch of the radial velocity 
curve. The second effect is a systematic difference between 
the observed and the computed phase lag between the ra- 
dial velocity curves and the light curves. Indeed, the phase 
lag of the BEP is positive and equal to Atpi = 0.055, while 
the phase lags for all our models (Tab. [3} are systemati- 
cally smaller and negative. A significantly better agreement 
is obtained for models with higher helium content, i.e. with 
X = 0.70. We also note that in spite of the model phase lags 
being all negative, the phase of maximum expansion velocity 
always precedes the phase of maximum brightness, just as in 
the case of the BEP. The negative value of A<^i results from 
the highly nonlinear shape of the light and radial velocity 
curves. 

(vi) The model for the BEP, that reproduces both the 
light and radial velocity variation best (has the lowest value 
of d a n) is the model with convective parameters of set S3 and 
with chemical composition of population I stars (X = 0.70, 
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Z = 0.01). Its absolute luminosity is 33 Lq and effective tem- 
perature is « 6970 K. Among the models with chemical com- 
position corresponding to population II stars, the best two 
models (with convective parameters of sets S2 and S3) have 
the same metallicity [Z = 0.01), luminosity (L = 32 Lq), 
and nearly the same effective temperature (« 6910 K). All 
our best matching models have very similar P2 /Po ratios, in 
agreement with interpretation that the bump is caused by 
the 2:1 resonance between the fundamental mode and the 
second overtone. 



4 NON-LINEAR MODEL SURVEY FOR BEPS 

At first, the BEP was recognised as a classical RR Lyrae 
star - its light variation and radial velocity curves mimic 
that observed in RR Lyrae stars pulsating in the fundamen- 
tal mode. Only its presence in the eclipsing binary system 
allowed to determine its mass and to reveal its true nature. 
Its mass is only half the typical RR Lyrae mass, too low to 
ignite helium. It cannot be a classical RR Lyrae star, its a 
product of mass transfer during the bina ry evolution. The 
incide nce rate of such stars estimated by IPietrzvhski et al.l 
l|2012l ) is only 2 per cent and only a fraction may be discov- 
ered through eclipses. Majority of such binaries are expected 
to be non-eclipsing systems. Consequently, direct estimate of 
mass of the pulsating component will not be possible. Then, 
the shape of the light/radial velocity curve may mimic that 
of RRab star, as it is the case for the BEP, and lead to 
wrong interpretation. The goal of this section is to conduct 
a large model survey of the light and radial velocity curves 
of BEPs, to check how to distinguish these stars from clas- 
sical RRab stars and to check whether the case of the BEP 
- RR Lyrae impostor - is a rule or an exception. 

We consider model sequences with three different 
masses corresponding to one-third - one-half the typical 
RR Lyrae mass, 0.2 M , 0.26 M Q (mass of the BEP) and 
0.3 Mq, and with different absolute luminosities, 20 Lq, 
25 L , 30 Lq, 35 Lq and 40 Lq. Models are computed along 
horizontal stripes of constant luminosity with a step of 100 K 
in effective temperature. Chemical composition is the same 
in all our models, X = 0.76 and Z — 0.01. Full model sur- 
vey was computed adopting convective parameters of set SI. 
In addition for 0.26 Mq we have computed model sequences 
adopting convective parameters of sets S2 and S3. 

In Fig. [9] we present the HR diagrams for models with 
M = 0.20 M (top panel) and M = 0.30 M (bottom panel, 
adopting convective parameters of set SI. For M — 0.26 Mq 
the HR diagram is plotted in the top left panel in Fig. [3] 
For the highest considered mass, the first overtone instability 
strip appears at the lowest luminosities. The models cover 
a wide period range from 0.3 up to 1.5 days. We focus our 
attention on models with Po below one day, as in RRab 
stars. As discussed in the previous section the 2:1 resonance 
is expected to shape the light and radial velocity curves of 
BEPs. For all considered masses, its loci (thick dotted lines 
in Figs. [5] and [3} crosses the instability strip at fundamental 
mode periods between 0.6 and 1 days. 

An analysis of the light and radial velocity curves and of 
the corresponding Fourier decomposition parameters clearly 
show that indeed, the 2:1 resonance with the second over- 
tone causes the progression of the light and radial velocity 



curve's shape. This is illustrated based on model sequences 
with M — 0.26 Mq (and different luminosities) and adopt- 
ing convective parameters of set SI. The bump progression, 
akin to that observed in classical Cepheids is best visible 
for radial velocity curves. The amplitude and the lowest or- 
der Fourier decomposition parameters are plotted in Fig. QJj] 
for model sequences with different luminosities. The Fourier 
parameters are plotted versus the P2/P0 period ratio. The 
progression of the Fourier coefficients clearly correlates with 
P2/P0 ratio and resembles ver y well the results obtained for 
classi cal Cepheid models (e.g. iBuchler. Moskalik fc Kovacsl 
1990) and predictions base d on analysis of amplitu de equa- 
tions for the 2:1 resonance l|Kovacs fc Buchlerll 19891 ). A clear 
drop in the pulsation amplitude is observed, followed by the 
decrease and minimum of the R21 ratio at somewhat larger 
values of P2/Po- In the case of the Fourier phase, a bell 
shape is apparent. The ascending part of the bell follows a 
universal progression - ip2i is a function of P2/P0 only (for 
P2/P0 > 0.52). The value of P2/P0 at which departure from 
this relation is observed and the highest possible value of (^21 
depends strongly on the luminosity. The correlation of the 
Fourier parameters progression with the P2/P0 ratio is lost 
if the Fourier parameters are plotted versus the period - see 
the blue curves in the right panel of Fig. 1141 which represent 
the same model sequences as plotted in Fig. 1101 No doubt, 
the 2:1 resonance is crucial in shaping the progression of the 
radial velocity curves. 

In Fig. [11] (right panel) we plot the radial velocity curves 
along a sequence of L — 25 Lq . In contrast to Fig. now a 
wide range of period ratios is probed, which clearly reveals 
the bump progression. The bump appears first in the middle 
of the ascending branch of the velocity curve, at the phase 
of maximum radius (P2/P0 ~ 0.54). Then, as period in- 
creases and period ratio decreases, the bump moves towards 
the phase of maximum expansion velocity (P2/P0 ~ 0.51). 
For period ratios below w 0.5 the bump is present on the 
descending branch and disappears for P2/P0 < 0.47. 

The bump progression is also visible for the J-band light 
curves (left panel of Fig. Ill|) . although in this case the in- 
terpretation is more difficult, because of additional features 
present in the light variation. At the highest effective tem- 
perature the light curve is roughly triangular. A trace of 
bump appears on the descending branch, just after the max- 
imum brightness, which next, moves towards the ascending 
branch and deforms the brightness maximum. For effective 
temperatures below 6700 K (P2/P0 < 0.52) a pronounced 
bump is present on the ascending branch. It is present even 
for the lowest temperature model, and in general, is a fea- 
ture of our models located in the centre and close to the 
red edge of the instability strip, even with period ratios as 
low as P2/P0 ~ 0.44 (not shown in Fig. Ill[) . although it 
is much less pronounced in such models. It is likely that 
other non-resonant effects play a role in shaping this bump. 
Another feature of the light curves, well visible in Fig. [TTJ 
is a nearly flat and long phase of high brightness for cooler 
models (T e g ^ 6200 K). At the lowest effective temperatures, 
two brightness maxima at similar level are present. Origin of 
such shape is most likely non-resonant and connected to the 
development of extensive zone of efficient convective transfer 
at the phase of maximum compression, which is sustained 
then for nearly half the pulsation cycle. 

Analysis of Figs. 1101 and ITTI indicates that the 2:1 res- 
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Figure 9. The HR diagrams for 0.20 M (top panel) and 0.30 Mq 
(bottom panel) models of set SI and with X = 0.76 and Z = 0.01. 
Long-dashed lines are lines of constant fundamental mode period 
(0.3d, 0.6d, l.Od and 1.5d). Loci of the 2:1 resonance with the 
second overtone and of the two half-integer resonances (Pq/Pi = 
1.5 and P0/P5 = 3.5) are also plotted. Triangles show the location 
of period-doubled models. 



onance influences the shape of the light and radial velocity 
curves in a wide range of period ratios, which we estimate to 
0.47 < P2/-P0 < 0.54. The border values of P2/P0 are plot- 
ted with thin dotted lines in Fig. [9] and labelled with P2/P0 
value. Clearly the 2:1 resonance is influential in a significant 
part of the instability strip and hence, is the most important 
effect shaping the light and radial velocity curves of BEPs. 

Another resonant effect which we find in our models 
is a period-doubling effect. For some models, we observe 
alternating deep and shallow minima/maxima in the light 
and radial velocity variation - see Fig. [12] for exemplary 
light curves. Consequently, pulsation repeats after two fun- 
damental mode periods. The behaviour is found only for 
some model sequences and is restricted to rather narrow pa- 
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Figure 10. Peak-to-peak amplitude and the lowest order Fourier 
decomposition parameters for the radial velocity curves for 
0.26 Mq models (set SI) of different luminosities (different sym- 
bols) plotted versus the linear P2/P0 period ratio. Pulsation pe- 
riod increases to the right. 



rameter ranges. For models adopting convective parameters 
of set SI, period doubling was found only in the least mas- 
sive, most luminous models. Their location is plotted with 

triangles in the top pa nel of Fig. [9] 

As analysed by iMoskalik fc Buchlerl (|l99fj ) period 
doubling in the hydrodynamic models is caused by the 
half-integer resonances. The effect is observed in sev- 
eral types of large amplitude pulsators including RV Tau 
variables (caused by the 5:2 resonance, P0/P2 = 2.5, 
IMoskalik fc Buchlerl Il99d ). RR Lyr stars showing the 
Blazhko effect (ca used by the 9:2 resonance, P0/P9 = 4.5, 
ISzabo et al]|2010h and i n BL Her stars (cau sed by the 3:2 
resonance, P0/P1 = 1.5. ISmolec et al.ll2012l ). In Fig. H] we 
plot the loci of two half-integer resonances that might be 
influential in the case of BEPs models (red lines). These 
are the 3:2 resonance with the first overtone (P0/P1 = 1.5) 
and the 7:2 resonance with the fifth overtone (P0/P5 = 3.5). 
To cause the period doubling, t he resonant overtone mod e 
cannot be damped too strongly (Moskalik & Buchlcr 1990). 
This is the case for both resonances. The fifth overtone is a 
trapped envelope mode i n the discussed models, fo r which 
the growth rate is high (|Buchler fc K ollath 200^). With- 
out the relaxation technique (e.g. IStellingwerJ Il974j ) how- 
ever, it is hard to identify which resonance is operational 
in these models. As the models are limited to narrow pa- 
rameter ranges we do not study this issue in more detail 
here. 
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Figure 11. /-band light curves (left panel) and radial velocity curves (right panel) for hydrodynamic models of 0.26 Mq computed along 
a line of constant luminosity, L = 25 Lq. Models adopt convective parameters of set SI. Effective temperature is given on the right 
hand side of the /-band curves, fundamental mode period and linear P2/P0 period ratio is given on the right hand side of the velocity 
curves. For clarity, the consecutive light curves were shifted by 0.2mag and radial velocity curves by 15kms _1 . In the latter case, short 
horizontal line segments mark the zero velocity level. Velocity curves were scaled by constant projection factor, p = 1.38. 



We now turn to the discussion of the overall properties 
of the model light and radial velocity curves for BEPs, in 
terms of their Fourier decomposition parameters. To pro- 
vide a feeling how reliable our non-linear code is in predict- 
ing the Fourier parameters, and to provide a background 
for discussing the parameters of BEPs, we have conducted a 
model survey for classical RR Lyrae stars, briefly mentioned 
in Section[2] results of which were presented in Fig. [2] These 
models adopt the convective parameters of set SI, but re- 
sults for other sets of convective parameters are very similar. 
Amplitude and -R21 amplitude ratio are modelled satisfacto- 
rily, a somewhat too low values of R31 are predicted. Fourier 
phases agree well with the observations, indicating that ma- 
jority of RR Lyrae stars have absolute luminosity around 
50 L . 

Fourier parameters from similar survey of BEP models 
are plotted in Figs. [13] (/-band light curves) and [14] (radial 
velocity curves) for all considered model sequences. Left pan- 
els in these Figures illustrate the effects of different mass of 
the models (with convective parameters of the models fixed 
to that of set SI), while right panels illustrate the effects 
of using different convective parameters in the modelling 
(sets SI, S2, and S3), while mass of the models is fixed to 
0.26 Mq. For each case, five model sequences with different 
absolute luminosities were computed. Successive models in 
each sequence were computed with a step of 100 K in effec- 



tive temperature. The first, shortest period model in each 
sequence is the hottest model located at the direct proxim- 
ity of the blue edge of the instability strip. It is clear that 
with increasing luminosity, model sequences shift towards 
longer periods. In all the Figures we also mark the location 
of the BEP (diamond) and for the /-band light curves we 
also plot the location of RRab and RRc (first overtone pul- 
sators) stars from the Galactic bulge for a reference. Such 
comparison is not possible in the case of radial velocity pa- 
rameters. Only for a few RR Lyrae stars coverage of radial 
velocity curves is good enough to allow determination of 
their Fourier parameters. Consequently, a meaningful com- 
parison is not possible at the moment. 

We first discuss the Fourier parameters for the /-band 
light curves displayed in Fig. [TS] As discussed before, we 
expect that the 2:1 resonance will have a significant effect 
in shaping the light curves. We first note that progressions 
of Fourier parameters with pulsation period are very sim- 
ilar for different model sequences displayed in the Figure. 
For example the R21 progressions have a wavy shape with 
local minima in the middle and at the edges of the progres- 
sion. R21 drops for the shortest and the longest period as a 
result of amplitude decrease at the edges of the instability 
strip. The decrease of R21 in the middle of the progression 
is a resonant effect (see Fig. I10[) . As resonance centre falls 
at different periods for different model sequences, the min- 
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Figure 12. Exemplary /-band light curves for models with 
period-doubling effect. 



imum of R21 also shifts in period. For example, for mod- 
els with M = 0.20 Mq the resonance centre is located at 
Po ~ 0.6 days for 20 Lq and at Po ~ lday for 40 L0 (top 
panel in Fig. [9} • 

The predicted peak-to-peak amplitudes of BEPs are be- 
low 0.6 mag in the /-band, which is smaller than the typ- 
ical amplitudes of RRab stars with periods shorter than 
~ 0.55 days. In this period range also the predicted am- 
plitude ratios are smaller than in RRab stars. We note how- 
ever, that the predictions for amplitudes and amplitude ra- 
tios should be treated with caution. Amplitudes and conse- 
quently amplitude ratios strongly depend on the convective 
parameters adopted in the model, which is well visible if 
one compares the results for the models of sets SI and S3 in 
the right panel of Fig. 1131 The lower the eddy viscosity (set 
S3), the higher become the amplitudes and amplitude ra- 
tios. Convective parameters should be calibrated to match 
the observational constraints, in particular the amplitude- 
period relation for a given group of pulsators. With only 
one BEP known at the moment the calibration is neces- 
sarily very rough. Fortunately, the Fourier phases do not 
depend strongly on the pulsation amplitude, which is clear 
from the analysis of the bottom right panels of Figs. 1131 and 
1141 The Fourier phases depend only weakly on the adopted 
convective parameters as well. Thus, the Fourier phases are 
the best and most reliable predictions of our modelling. 

The progressions for the Fourier phases (bottom panels 
of Fig. 1 13f) are very interesting. Note that for ^31 the model 
sequences were plotted twice, shifted by 2-7T. The 2-7T ambigu- 
ity in the values of the Fourier phases should be kept in mind 
while analysing the Figures. We note that 1^21 progressions 
strongly depend on the mass and on the absolute luminosity 
adopted in the models. Thus, </P2i may provide useful con- 
straints on the physical parameters of BEPs discovered in 
the future. Although for the majority of our models </>2i sig- 
nificantly differs from the RRab progression, values typical 
for RRab stars are also possible, particularly for models lo- 
cated on the hot side of the instability strip. Consequently 
if 21 is not the best indicator of BEP. It is the 1/P31 Fourier 
phase that may be used to distinguish BEPs from RR Lyrae 
stars. Clearly, the values of 9331 predicted for BEPs differ sig- 
nificantly form the values typical for RR Lyrae stars. Only 
very close to the blue edge of the instability strip (and this 
is the case for the BEP) the values of 1^31 are similar for the 



two groups. Through the majority of the instability strip 
<£3i is larger in the BEPs models. For cooler models the !/?3i 
value increases. In the middle of the instability strip the dif- 
ference between ipz\ as predicted for BEPs and the RRab 
progression is around n rad. For lower masses the difference 
becomes smaller, but even for 0.2 Mq models the difference 
is about 2 rad. 

We also note that at shorter periods BEP models are 
also different from RRc stars. The predicted amplitudes are 
larger than those of RRc stars, and the same is true for the 
-R21 ratio. The modelled Fourier phases, both <p2i and 9531, 
fall in between the RRab and RRc progressions. The phases 
are larger than in RRab stars, but smaller than in RRc stars. 

In the case of the radial velocity curves, the Fourier 
parameters displayed in Fig. [14] should be treated as pre- 
dictions, to be compared with future observations of BEPs. 
It is hard to pinpoint the expected differences between the 
Fourier parameters for BEPs and those for RRab stars, as 
radial velocity observations of RRab stars are currently very 
scarce. The 1^21 values seem to be the most interesting ones 
for radial velocity curves. They depend only slightly on the 
adopted convective parameters. The bell shape is a result of 
the bump progression, and this is the most important differ- 
ence between BEPs and RRab stars. No bump progression 
is expected in the latter case. The presence of a bump on a 
radial velocity curve is thus a good indicator for a BEP. The 
situation is more difficult in the case of the light curves, as 
additional features (bumps, shoulders) which may result for 
example from the shock pr opagating in the atmosphere (e.g. 
iBono fc Stellingwerj|l994r ). are commonly present in RRab 
stars. We also note that the <p2i values for the radial velocity 
curves might be used as an indicator of mass and absolute 
luminosity. 



5 DISCUSSIONS AND CONCLUSIONS 

In Section [3] we have conducted a model survey to repro- 
duce the observed properties of the BEP. Our best model 
adopts the convective parameters of set S3, and a chemical 
composition typical for population I stars, i.e. X = 0.70 and 
Z = 0.01. The model is located at the blue edge of the in- 
stability strip, its absolute luminosity is 33 Lq and effective 
temperature is T e g s=s 6970 K. The model light curve fits the 
observed light variation very well. Good agreement was also 
obtained for the radial velocity curve (Fig. [HJ . We have re- 
produced the bump on the ascending branch of the radial 
velocity curve of the BEP, and traced its origin to the 2:1 
resonance between the fundamental mode and the second 
overtone, P2/P0 = 0.5. 

The models with a chemical composition typical for 
population II stars (X = 0.76) also reproduce the observa- 
tions of the BEP well and have nearly the same luminosity 
and effective temperature. We note that although the mass 
of the BEP is very low and its pulsation resembles that of 
an RR Lyrae star, the BEP is the product of mass transfer 
in the binary system. 5.4 Gyr ago t he progenitor of the BEP 
was a 1.4 M© star, as models of IPietrzvnski et alj (|2012| ) 
indicate and thus, the chemical composition typical for pop- 
ulation I stars is very likely the correct one. 

We also considered models with other values of convec- 
tive parameters. The resulting properties of the best match- 
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Figure 13. Fourier decomposition parameters for the /-band light curves of non-linear models of BEPs. Left panel: models of different 
masses, 0.20 Mq, 0.26 Mq and 0.30 Mq along five horizontal stripes with different luminosity (20 — 40 Lq, different symbols). Right 
panel: Models of 0.26 Mq and different luminosity and adopting different values of convective parameters (sets SI, S2 and S3). Diamond 
indicates the position of the BEP. For comparison, data for classical RR Lyrae stars from the Galactic bulg e from the OGLE catalogue 
are plotted (grey dots for RRab stars, magenta dots for RRc stars, each third plotted; ISoszvnski et al]|2011f ). Sine convention is used for 
Fourier phases. 



ing models for the BEP are then very similar to those quoted 
above (see Tab. [3]) and therefore depend only slightly on the 
detailed treatment of convection in the model. 

We note that the phase lag between the radial veloc- 
ity curve and the light curve is systematically smaller in 
our models than that observed in the BEP. More computa- 
tions are needed to find out how this discrepancy might be 
removed. The current models indicate that a slight helium 
enrichment might be needed. Also effects of turbulent con- 
vection, not considered in this study (e.g. turbulent pressure, 



overshooting) need to be checked. A more detailed modelling 
of the BEP is planned in the near future. 

The effective temp erature of the BEP estimated by 
IPietrzvriski et al.1 (|2012l ) is « 350 K larger (T cff = 7320 ± 
160 K) than in our best matching model. Consequently, also 
the absolute luminosity is higher (L = 46.2 ± 9.3 L©). The 
estimates of effective temperature and absolute luminosity 
of the BEP not only differ from our results, but are also in 
conflict with stellar pulsation theory. They place the BEP 
outside the instability strip, as the analysis of Fig. [3] clearly 
indicates. This conclusion is independent of the convective 
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Figure 14. Fourier decomposition parameters for the radial velocity curves of non-linear models of BEPs. Left panel: models of different 
masses, 0.20 Mg ,0.26 Mg and 0.30 Mg along five stripes of different luminosity (20 — 40 L0). Right panel: Models of 0.26 Mg and different 
luminosity and adopting different values of convective parameters (sets SI, S2 and S3). Diamond indicates the position of the BEP. Sine 
convention is used for Fourier phases. 



parameters and chemical composition adopted in the model 
calculations. 

The estimate of the e ffective temperature of the BEP 
bv lPietrzvriski et alj l|2012l ) was based on the assumed T2 = 
5000 K effective temperature of the secondary component. 
It is a typical temperature for a 1.67 Mg giant, but it is not 
an actual measurement. Then, the effective temperature of 
the pulsating component of the binary was calculated us- 
ing a temperature ratio (Ti/Tb = 1.464 ± 0.010) obtained 
from the analysis of the light curve of the binary. To check 
how the assumption on T2 affects the resul ts we have re- 
done the analysis of iPietrzvhski et al.l (|2012l ) but adopting 
lower values of Ti. It turned out that a decrease of T2 by 



100 K leads to a decrease of T\ by w 240 K. Thus, only a 
slight decrease of the effective temperature of the secondary 
component, by « 150 K, brings the effective temperature of 
the BEP into agreement with the effective temperature re- 
sulting from our best models. The absolute luminosity is still 
somewhat too high (L ~ 38 Lq) as compared to our models, 
which is caused by slightly different estimates of the radius. 

The radius of the BEP estimated by IPietrzvhski et all 

l|2012D is R = 4.24 ± 0.24 R , and, as we checked, is inde- 
pendent on the assumed value of T2. The values inferred 
from our models (see Fig. 0and Tab. [3} are lower by up to 
1.5(7 (1.5 a for models with Z — 0.01 and 1.2 a for mod- 
els with Z — 0.001). Multicolour photometry of the binary 
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system may help to determine the radius of the BEP more 
accurately and this way tell us how severe the discrepancy 
is. 

We have conducted a large non-linear model survey for 
the new class of variable stars. The aim was to provide the 
criteria to distinguish BEPs from classical RR Lyrae pul- 
sators. The resonant bump progression caused by the 2:1 
resonance between the fundamental mode and the second 
overtone is the main factor expected to determine the shapes 
of the light and radial velocity curves of BEPs. This is not 
the case for classical RR Lyrae stars. In the latter case, the 
2:1 resonance falls close to the red edge of the instability 
strip and hence, no bump progression is observed in clas- 
sical RR Lyrae stars. Thus, the presence of a bump may 
indicate the binary evolution of the RR Lyrae impostor. We 
note that the location of the bump on the radial velocity 
curve is a measure of the P2/P0 ratio, which provides an 
excellent constraint on the stellar models. In the case of the 
light curves, one may expect pronounced bumps on the as- 
cending branch and long and flat, or long and double-peaked 
phases of high brightness (Fig. Ill[) . In some of our models 
with the lowest M/L ratio we find a period-doubling effect 

(Fig.rju. 

The other excellent discriminant of the new class of pul- 
sating stars is the value of the ipzi Fourier phase for the I- 
band light curve. Except the models located extremely close 
to the blue edge of the instability strip, the expected values 
of 1^31 differ significantly from the typical values for RRab 
stars. The values of ip^i on the other hand, might be use- 
ful indicators of stellar mass and absolute luminosity. We 
also note that BEPs are expected to chan ge (decrease) their 
periods very fast l|Pietrzvnski et al"1l2012r ). 

Our models clearly indicate that the BEP with its light 
and radial velocity curves strongly mimicking the typical 
RR Lyrae curves, might be a rare case. For the majority of 
BEPs, significant departures from the observational proper- 
ties of RR Lyrae stars are expected. 



on the psk computer cluster in the Copernicus Astronomical 
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